(N 






% 



UK/97-23, ADP-97-55/T281 

Dec. 1997 
hep-ph/9712483 



Lattice Calculation of the Strangeness Magnetic 
Moment of the Nucleon 

"S" 

On 

5} ■ S.J. Dong a , K.F. Liu a , and A.G. Williams'' 

O ! a Dept. of Physics and Astronomy, Univ. of Kentucky, Lexington, KY 40506 

Q 



b Special Research Cneter for the Subatomic Structure of Matter and Department of 



Physics and Mathematical Physics, University of Adelaide, Australia 5005 



Abstract 



> 

m 

oo 

We report on a lattice QCD calculation of the strangeness magnetic moment 
of the nucleon. Our result is G S M (0) = —0.36 ± 0.20. The sea contributions 

^. , from the u and d quarks are about 80% larger. However, they cancel to a large 

extent due to their electric charges, resulting in a smaller net sea contribution 

r S^ • of —0.097 ± 0.037^tv to the nucleon magnetic moment. As far as the neutron 

to proton magnetic moment ratio is concerned, this sea contribution tends to 
cancel out the cloud-quark effect from the Z-graphs and result in a ratio of 
—0.68 ± 0.04 which is close to the SU(6) relation and the experiment. The 
strangeness Sachs electric mean-square radius (rf)E is found to be small and 
negative and the total sea contributes substantially to the neutron electric form 

/\ ' factor. 



PACS numbers: 12.38.Gc, 13.40.Fn, 14.20.Dh 



The strangeness content of the nucleon has been a topic of considerable recent 
interest for a variety of reasons. The studies of nucleon spin structure functions in 
polarized deep inelastic scattering experiments at CERN and SLAC M, combined 
with neutron and hyperon /3 decays, have turned up a surprisingly large and negative 
polarization from the strange quark. In addition, there is a well-known long-standing 
discrepancy between the pion-nucleon sigma term extracted from the low energy pion- 
nucleon scattering ||] and that from the octect baryon masses ||. This discrepancy 
can be reconciled if a significant ss content in the nucleon 0, |3[] is admitted. 

This naturally leads to the question that if the strange quarks can contribute 
substantially in the axial-vector and scalar current matrix elements, how important 
are they in other matrix elements involving the vector, pseudoscalar, and tensor 
currents. The case of the vector current matrix element (iV|s7 M s|./V) is especially 
interesting. If the strange magnetic moment (m. m.) is large, it is likely to spoil 
the nice SU(6) prediction of the neutron to proton m. m. ratio of —3/5 which lends 
credence to the valence quark picture. On the other hand, if it is small one would 
like to understand why it should be different from the axial-vector and scalar cases. 

To address some of these issues, an experiment to measure the neutral weak mag- 
netic form factor Gfj via elastic parity-violating electron scattering at backward an- 
gles was recently carried out by the SAMPLE collaboration ||. The strangeness 
magnetic form factor is then obtained by subtracting out the nucleon magnetic form 
factors G P M and G n M . The reported value is G S M (Q 2 = O.lGeV 2 ) = +0.23 ± 0.37 ± 
0.15 ± 0.19, where the last error is due to an uncertainty associated with axial radia- 
tive corrections. Future experiments have the promise of tightening the errors and 
isolate the radiative corrections so that we can hope to have a well-determined value 
and sign for G S M (0). 

Theoretical predictions of G S M (0) vary widely. The values from various models 
and analyses vary from —0.75 ± 0.30 in a QCD equalities analysis to +0.37 in 
an SU(3) chiral bag modelf7|. While a few give positive values 0, ||, most model 
predictions are negative with a typical range of —0.25 to —0.45 || £| [K], [11], [12], [13], |H 



Summaries of these predictions can be found in Refs. p, |i4}|. A similar situation exists 
for the strangeness electric mean-square radius (t 2 )e- A number of the predictions 



are positive || |i~3f , while the others are negative || [II], |TT| [12], |H}| . Elastic e p and 



e He parity- violation experiments are currently planned at TJNAF [I5| to measure 
the asymmetry Am at forward angles to extract the strangeness electric mean-square 
radius. Hopefully, they will settle the issue of the sign of (v 2 )e- 

In view of the large spread of theoretical predictions for both G S M (0) and (t 2 )e 
and in view of the fact that the experimental errors on G S M (0) are still large, it is 
clearly important to perform a lattice calculation of this quantity since this is a first 
principles theoretical approach. We present our calculation from lattice QCD in the 
hope that this will shed some light on these quantities. Our previous results on flavor- 
singlet quantities which involve the so-called "disconnected insertions" (DI) for the 
sea quarks in addition to the "connected insertions" (CI) for the valence and cloud 



quarks |L6], [L7] reveal that the sea quark contribution to the flavor-singlet g A from 



the DI is negative and the magnitude large enough (e.g. the strangeness polarization 
As = 0.12±0.01) to cancel the positive CI contribution to a large extent. This results 
in a small g\ at 0.25 ± 0.12, which is in agreement with the experimental results [pi. 



Similarly, the calculated ratio y = (N\ss\N) / (N\uu + dd\N) = 0.36 ± 0.03 Jl7j gives 
the right amount of strangeness content to resolve the ttNu puzzle we alluded to 
earlier. Given these reasonably successful estimates of strangeness in the axial- vector 
and scalar channels, we feel that it should yield meaningful results in the vector 
current as well. In particular, we would like to understand why the SU(6) valence 
quark picture fails badly in the flavor-singlet axial-vector and scalar cases and yet 
gives an apparently good prediction in the neutron to proton m. m. ratio - a yet 
unresolved puzzle in low-energy hadron physics. 

The lattice formulation of the electromagnetic form factors has been given in de- 
tail in the past |18 , [H|. Here, we shall concentrate on the DI contribution, where 



the strangeness current contributes. In the Euclidean formulation, the Sachs electro- 
magnetic form factors can be obtained by the combination of two- and three-point 
functions 

G% a N (t,p) =J2^ lf *( Q \x a (x)x a (0)\0) (1) 

X 

G a N % N (t f ,p,t,q) = £ e-^f + ^(0\ X a (x f )V,(x)x P (0)\0), (2) 

x f ,x 

where x a is the nucleon interpolating field and V^(x) the vector current. With large 
Euclidean time separation, i.e. tf — t » a and t » a, where a is the lattice spacing, 

rf a Gf VjN (t f , 0, t, q) G%%(t, 0) e ijk q k 2 



Gffif(t f ,0) G% a N {t,q) E q + m 

^E a G%y 4N (tf, 0, t, q) G^ N (t, 0) 
G a N a N (tf,0) G a N a N (t,q) 




G E (q 2 ), (4) 



where 1^ = I Q * 1 , and T E = 

We shall use the conserved current from the Wilson action which, being point- 
split, yields slight variations on the above forms and these are given in Ref. [18, 19 . 
Our 50 quenched gauge configurations were generated on a 16 3 x 24 lattice at (3 — 6.0. 
In the time direction, Dirichlet boundary conditions were imposed on the quarks to 
provide larger time separations than available with periodic boundary conditions. We 
also averaged over the directions of equivalent lattice momenta in each configuration; 
this has the desirable effect of reducing error bars. Numerical details of this procedure 



are given in Refs. [JT9], [2l[] . The dimensionless nucleon masses M^a for k = 0.154, 
0.152, and 0.148 are 0.738(16), 0.882(12), and 1.15(1) respectively. The corresponding 
dimensionless pion masses m^a are 0.376(6), 0.486(5), and 0.679(4). Extrapolating 
the nucleon and pion masses to the chiral limit we determine k c = 0.1567(1) and the 
dimensionless nucleon mass at the chiral limit to be 0.547(14). Using the nucleon mass 



to set the scale, which we believe to be appropriate for studying nucleon properties p0, 



|2~T| , [16], [T7||, the lattice spacing a 1 = 1.72(4) GeV is determined. The three k's then 
correspond to quark masses of about 120, 200, and 360 MeV respectively. 

The strangeness current s'-f^s contribution appears in the DI only. In this case, 
we sum up the current insertion time slice t from the nucleon source to the sink in 
Eqs.((|) and (||) to gain statistics flB], |I7|. This leads to const +tfGE,di S (q 2 ) for Eq. ([|). 



For Eq. (|3]), we average over the three spatial components s'fiS and obtain const + 



tf e + m ^ M ' dis (ff 2 )- Similar to our previous studies of As [[nj and (iV|ss|iV) fI7fl , we 



use 300 complex Z 2 noises [Q and 100 gauge configurations to calculate the sea quark 
contribution (DI) with k = 0.148, 0.152 and 0.154. In calculating the strange current, 
we have considered the correlation between the quark loop with k s = 0.154 and the 
valence quarks at k v = 0.148, 0.152, and 0.154. The ratio in Eq. (|3|) with the sum in 
t and average in Vi, which leads to the expression const + tfGM,dis(q 2 ), is plotted in 
Fig. 1 as a function of tf for |g| = 2tc /La. Then GM,dis(q 2 ) from the DI is obtained 
from fitting the slopes in the region £/ > 8 where the nucleon is isolated from its 
excited states with the correlation among the time slices taken into account jT(|. The 
resultant straight-line fits covering the ranges of tf with the minimum \ 2 are plotted 
in Fig. 1. Finally, the errors on the fit, also shown in the figure, are obtained by 
jackknifing the procedure. To obtain the physical G s M (q 2 ), we extrapolate the valence 
quarks to the chiral limit while keeping the sea quark in the loop at the strange 
quark mass (i.e. k s = 0.154). It has been shown in the chiral perturbation theory 
with a kaon loop that G S M (0) is proportional to mx, the kaon mass [E3]. Thus, we 



extrapolate with the form C + D^Jrh + m s where m is the average u and d quark 
mass and m s the strange quark mass to reflect the mx dependence. This is the same 



form adopted for extracting (N\ss\N) in Ref. [|17 , which also involves a kaon loop in 
the chiral perturbation theory. 
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Figure 1: The ratio of Eq. (|3|) as a function of tf so that the slope is Gm, dis(<? 2 ) at 
q — 2n/La. The sea quark is fixed at k s = 0.154, the strange quark mass, and the 
valence quark masses are at k v = 0.148, 0.152, and 0.154. M is the fitted slope. 

Plotted in Fig. 2(a) is the extrapolated G s M (q 2 ) at 4 nonzero q 2 values. The errors 
are again obtained by jackknifing the extrapolation procedure with the covariance 



matrix used to include the correlation among the three valence «'s. In view of the 
fact that the scalar current exhibits a very soft form factor for the sea quark (i.e. 
9s,dis(c[ 2 )) which has been fitted well with a monopole form [17], we shall similarly 
use a monopole form to extrapolate G s M (q 2 ) with nonzero q 2 to G^(0). Indicated 



as o in Fig. 2(a), we find G S M (0) = —0.36 ± 0.20. Again, the correlation among the 
4 q 2 are taken into account and the error is from jackknifing the fitting procedure. 
This is consistent with the recent experimental value within errors (see Table 1). The 
monopole mass is found to be 0.58 ± 0.16 of m^. To explore the uncertainty of the q 2 
dependence, we also fitted G s M (q 2 ) with a dipole form and found G S M (0) = — 0.27±0.12 
with a dipole mass m^/m-Tv = 1.19 ± 0.22. Similar results are obtained for u and d 
quarks with monopole fits. They turn out to be G^ dis (0) = —0.65 ± 0.30, which is 
about 1.8 times the size of G S M (0). These are tabulated in Table 1. 
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Figure 2: (a) Strange magnetic form factor G s M (q 2 ). G S M (0), indicated by o, is ob- 
tained from a monopole fit. (b) Strangeness electric form factor G s E (q 2 ). The solid 
line is a fit with the monopole form shown in the figure and the dashed line is obtained 
with the monopole mass thm from G s M (q 2 ) in (a). 



Now, we are ready to address the question of why the SU(6) relation is badly 
broken in the scalar current (e.g. F s , D s ) and axial current (e.g. g^) and yet is so 
good for the neutron-proton magnetic moment ratio /x n //x p . The lattice calculations 



for the scalar [|T7J and axial [HJ currents reveal the fact that the SU(6) breaking comes 
from both the sea quarks in the DI and the cloud quarks in the CI. We shall see how 
these degrees of freedom play out in the case of the m. m. We first plot in Fig. 3 
the ratio fi n /^ P for the CI part (shown as o) as a function of the valence quark mass. 
We see that when the quark mass is near the charm region {m q a at 0.55 corresponds 
to m q ~ 1 GeV), the ratio is close to the SU(6) prediction of -2/3. This is quite 
reasonable as this is in the non-relativistic regime where one expects SU(6) to work 
well. As the quark mass comes down to the strange region [m q a = 0.07), the ratio 
becomes less negative. Extrapolated to the chiral limit, the ratio is —0.616 ± 0.022 
which deviates from the SU(6) prediction by 8 %. We understand this deviation as 



mainly due to the cloud quark effect in the Z-graphs. As we switch off these Z-graphs 



in a valence approximation [|T6|, |l7j], the ratio (plotted as o in Fig. 3) becomes closer 
to the SU(6) value which resembles the non-relativistic case. Similar behaviors were 
observed for the scalar and axial matrix elements [[17], |16| . Now we add the sea quark 
contribution from the DI to give /jL dis = (2/3G^ dis (0) - l/3G d Mdis (0) - l/3G s M (0))fi N 
to the CI and find that it tends to cancel the cloud effect and bring the ratio back 
to be similar to what the valence approximation predicts. For the G S M (0) at various 
K v , we use the G s M (0)/G^. Idis (0) ratio from the chiral limit to obtain it from the 
^Mdis(O) a t each k v . At the chiral limit, when the total sea contribution fidis = 
— 0.097±0.037yU7v is added to the CI, the fi n /^ P ratio then comes down to — 0.68±0.04 
which is consistent with the experimental value of 0.685. We note that the fi n /fi p 
ratio for the full result (•) is more negative at the chiral limit compared with those 
at other m q a. This has to do with the fact that the CI employs the linear quark mass 
extrapolation, as do other observables for the CI [IIB, O, |21[ , whereas the DI uses the 



m„ dependence for the chiral extrapolation as mentioned above. From this analysis, 



we see that although the individual G^ dis (0), G d M dis (0), and G S M (0) are large, their 
net contribution /i dis = —0.097 ± 0.037/iAr to // n and /j, p is much smaller because of 
the partial cancellation due to the quark charges of u, d, and s. The sea contribution 
turns out to be further canceled by the cloud effect to bring the fi n /n P ratio close 
to the experimental value and the SU(6) relation. Barring any known symmetry 
principle yet to surface, this cancellation is probably accidental and in stark contrast 
with the nNcr term and flavor-singlet g° A where the cloud and sea effects add up to 
enhance the SU(6) breaking [17, [16 . 
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Figure 3: Neutron to proton m. m. ratio fi n /fi p as a function of the dimensionless 
quark mass m q a. The o indicates the valence result. The o is the result for the CI 
and • indicates the full result with the inclusion of the sea from DI. 



The sea contribution from the u, d, and s quarks G^' Idis (q 2 ) and G s M (q 2 ) are 
added to the valence and cloud part in the CI, GM,con(<? 2 ), to give the full G p M (q 2 ) 
and G%[(q 2 ). They are plotted in Fig. 4(a) and 4(b) and indicated by •. Also plotted 
are the GM,con(q 2 ) (denoted by o) and the experimental fits (in solid line). We see 
from Fig. 4 and Table 1 that \x v and //„ are smaller than the experimental results 
by ~ 6% in absolute values. This is presumably due to the systematic errors of the 
finite volume, finite lattice spacing, and the quenched approximation. We should 
point out that in the earlier discussion of the neuton to proton m. m. ratio /x n /// p , 
the systematic errors are expected to be cancelled out in the ratio to a large extent. 
Our conclusion of the ratio \i n j \i v in the preceding paragraph is thus based on this 
assumption. 
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Figure 4: (a) Proton magnetic form factor G p M (q 2 ). o indicates the result from the 
CI. They are shifted slighted to the left in — q 2 to avoid overlap with the full result 
which is shown as •. The solid line is the fit to the experiment [24]. (b) the same as 
in (a) for the neutron form factor G r l I {q 2 ). 



A similar analysis is done for the strange Sachs electric form factor G s E (q 2 ). This 
is plotted in Fig. 2 (b). We see that G S E (0) is consistent with zero as it should be. 
This serves as a test of the stochastic noise estimation with the Z 2 noise. We fitted 
G s E (q 2 ) with the form G s E (q 2 ) = / 4-/(1 - q 2 /m 2 M ) (solid line in Fig. 2(b)). The 



resultant electric mean-square radius (r 2 )j 
This is shown in Table 1. 
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dG%(q 2 ) 
dq 2 



-0.061 ± 0.003 fm 2 



In view of the large errors, we also plot the above form for G s E (q 2 ) with the 
monopole mass thm taken from G s M (q 2 ) which is shown by the dashed line. This 
gives (r 2 s ) E = -0.16 ± 0.06 fm 2 with x 2 /N DF = 0.24. This shows that that the 
uncertainty in the fitting can be as large as a factor of two. Nevertheless, (r 2 ) E is 
relatively small. This small negative value in (r 2 ) E and large negative G S M (0) are 
consistent with the kaon loop picture |11] and VMD [12| but is inconsistent with 



most of the other model predictions [0, WA 



Since the DI of u and d quarks are slighter larger than that of the s quark, the 
total sea contribution G E:dis (q 2 ) = 2/3G u Edis (q 2 ) - l/3G d Edis (q 2 ) - l/3G s E (q 2 ) adds a 
small positive value to the valence and cloud part G EtCon (q 2 ) in the CI. The proton 
G p E (q 2 ) and neutron G E (q 2 ) are plotted in Figs. 5(a) and 5(b) respectively. We see 
that the CI part G P E con (q 2 ) (shown as o in Fig. 5(a)) gives the main contribution 
in proton. G E ^ s (q 2 ) adds only a little change to it. The resultant dipole fit gives a 
dipole mass of 0.857 ±0.031 GeV (Table 1). This is consistent with the experimental 
dipole mass of 0.842 GeV. In the case of the neutron, since G Econ (q 2 ) |19[ itself (o 
in Fig. 5(b)) is small, the sea contribution G Et dis(q 2 ) becomes a sizable part of the 
total G E (q 2 ) (• in Fig. 5(b)). We see that when the sea is included we have a 
reasonably good match with the experimental results (solid line in Fig. 5(b)). The 
total mean square charge radius of —0.123 ± 0.019fm 2 is obtained from fitting with 
the form G E (q 2 ) = /t4-/(1 - q 2 /M 2 D f/{l - h.Qq 2 /AM 2 N ) which has been used to fit 



the experimental results |24 
of -0.127fm 2 . 



This is consistent with the experimentally fitted result 
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Figure 5: (a) Proton electric form factor G p E (q 2 ). o shows the result from the CI. 
They are shifted slighted to the left in — q 2 to avoid overlap with the full result which 
is shown as •. The solid line is the fit to the experiment [24]. (b) the same as (a) for 
the neutron form factor G E (q 2 ). 

In summary, we have calculated the s and u, d contributions to the electric and 
magnetic form factors of the nucleon. The individual m. m. and electric form factors 
from the different flavors in the sea are not small, however there are large cancellations 
among themselves due to the electric charges of the u, d, and s quarks. We find that a 
negative G S M (0) leads to a total negative sea contribution to the nucleon m. m. which 
cancels the cloud effect to make the fi n / 'fi p ratio consistent with the experiment. 
We also find G s E (q 2 ) positive and leads to a postive total sea contribution to the 
neutron electric form factor G E (q 2 ). Future calculations are needed to investigate the 
systematic errors associated with the finite volume and lattice spacing as well as the 
quenched approximation. 
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Table 1: Strangeness and proton- neutron m. m. and charge radii in comparison with 
experiments. 

Lattice Experiments 

G S M (Q 2 = O.lGeV 2 ) = 0.23 ± 0.37 ± 0.15 ± 0.19 



2.79 fjt N 

-l.9lfl N 

-0.685 



GU0) 


-0.36 ±0.20 


^M,dis(0) 


-0.65 ±0.30 


fidis 


-0.097 ± 0.037^ 


th 


2.62 ± 0.07 [in 


fin 


-1.81 ± 0.07 n N 


fin/ fip 


-0.68 ±0.04 


{rl)E 


-0.061(3) 0.16(6) fm 2 


(r 2 ) P E 


0.636 ± 0.046 fm 2 


(r 2 ) n E 


-0.123 ± 0.019 fm 2 



0.659 fm 2 [|4 
-0.127 fm 2 |24| 
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